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Q . Abstract 
(N ■ 

Let a three-dimensional isotropic elastic body be described by the Lame system 
^ • with the body force of the form F{x,t) = ip{t)f{x), where if is known. We consider 

I the problem of determining the unknown spatial term f{x) of the body force when 

I the surface stress history is given as the overdetermination. This inverse problem is 

ill-posed. Using the interpolation method and truncated Fourier series, we construct a 
. regularized solution from approximate data and provide explicit error estimates. 
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1 Introduction 

m 

1> • Let Q = (0, 1) X (0, 1) X (0, 1) represent a three-dimensional isotropic elastic body and let 
T > be the length of the observation time. For each x := {xi,X2,X3) G fl, we denote by 
! u{x,t) = {ui{x,t),U2{x,t),U3{x,t)) the displacement, where uj is the displacement in the 
Xj-direction. 

As known, u satisfies the Lame system (see, e.g., [HE]) 



o 



— -/iA^i-(A + /i)V(div(u)) = F, xen,teiO,T), (1) 

^ ' where F{x, t) = {Fi{x, t), F2{x, t), F^^x, t)) is the body force and div(u) = V -m = dui/dxi + 
du2/dx2 + dus/dx^. The Lame constants A and /x satisfy ^ > and A + 2/i > 0. The system 
([1]) is associated with the initial condition 

(ui(x,0),M2(a;,0),U3(a;,0)) = {gi{x), g2ix), g^ix)), x E il, 

(a;, 0), ^(x, 0), ^(x, 0)^ = {hi{x), h2{x), hsix)), x e Q, ^ ^ 



dt ' ' dt ' ' dt 
and the Dirichlet boundary condition 

(mi(x, t),U2{x, t),U3{x, t)) = (0, 0, 0), xedQ,te (0, T), (3) 
namely the boundary of the elastic body is clamped. 
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The direct problem is to determine u from u{0,x),Ut{0,x) and F. We are, however, 
interested in the inverse problem of determining both of {u,F). Of course, to ensure the 
uniqueness of the solution we shall require some additional information (the overdetermina- 
tion). Similarly to [21 16], we shall assume that the surface stress is given on the boundary of 
the body, i.e.. 



where n = (rii, n2, n^) is the outward unit normal vector of dfl and the stresses a and r are 
defined by 



In 2005, Grasselli, Ikehata and Yamamoto p] showed that the body force of the form 
F{x,t) = ip{t)f{x) is uniquely determined from (IIM]) provided that (p G C^([0,T]) is given 
such that ^p{0) ^ and the time of observation T > is large enough. In spite of the 
uniqueness, the problem of determining the spatial term / is still ill-posed, i.e. a small error 
of data may cause a large error of solutions. Therefore, it is important in practice to find 
a regularization process, namely to construct an approximate solution using approximate 
data. 

Recently, the regularization problem was solved partially in [6J, where a regularized so- 
lution for the time-independent term / is produced using further information on the final 
condition u{x,T). The final condition plays an essential role in [S] since it enables the au- 
thors to find an explicit formula for the Fourier transform of /, and then use this information 
to recover /. 

It was left as an open problem in j6j (see their Conclusion) to find a regularization 
process without using this technical condition. The aim of the present paper is to solve 
this problem completely, i.e. to find a regularization process of / using only the data in 
(PH]). We follow the interpolation method introduced in [7] where the authors constructed a 
regularized solution for the heat source of a heat equation. More precisely, lacking the final 
condition, we are only able to find an approximation for the Fourier transform /(^) with |^| 
large. The idea is that because /(^) is an analytic function (since / has compact support), 
we can use some interpolation process to recover f{^) with |,^| small. 

The paper is organized as follows. In Section 2, we shall set some notations and state our 
main results. Then we prove the uniqueness in Section 3 and the regularization in Section 
4. Finally, in Section 5 we test our regularization process on an explicit numerical example. 




(4) 




2 Main results 



Recall that our aim is to recover the spatial term 



fi^) = {fl{x), f2{x), f3{x)),X G Q, 
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of the body force F{x,t) = (f{t)f{x) from the system ([Ull]). The Lame constants always 
satisfy /i > and A + 2/i > 0, and the data / = {ip, X, g, h) is allowed to be non-smooth, 

/ G (Li(0,T), {L\Q,T,L\dm\ {L\n))\ {L\n)f) . 

For e = (^1, 6, es), C = (Ci, C2, Cs) G C3, we set e ■ C = 6Ci + 6C2 + iCl = and 
l^lg = .^|. For a e and X G R^ denote 

(^(ajx) = Gl{a,x) = Gl{a,x) = Gl{a,x) 

= cos (aiXi) cos (023^2) cos (03X3) , 
(a, x) = (a, x) = — sin (aiXi) sin (02X2) cos (03X3) , 
Gl{a,x) = (a, x) = — sin (aiXi) cos (02X2) sin (03X3) , 

G3(a,x) = G2 (tt, x) = — cos (ttiXi) sin (0:2X2) sin (03X3) . 

Sometimes we shall write G;^ instead of G;^(a,x) if there is no confusion. 
We start with the following lemma. 

Lemma 1. Ifu G {G\[0,T]; L\n))r]L^{0,T; H\n))f and f e {L\n)f satisfy the system 
(Qlj^ with data I = {(f, X, g, h), then 



E,,{I) {a) + EUa) , E^.il) {a) + EU 



a] 



n 



Di[I){a) D2{I){a) 



for all a = (ai, 0:2, 03) G such that 

\a\l = -{al + al + al)>0 and (/) (a) ^ 0, D2 (/) {a) ^ 0, 

where 



sinh v/(A + 2//)|a|o(T-t) 

(/)(«) = -\a\lj^{t) ^ . 

cosh ( Y (A + 2/i) |q;|qT 







/• smh(y^|a|o(T-t) 

D2 (/)(«) = -\a\l ifit) ^ ^dt, 

{ cosh(^^|a|oTj 

Eij{I){a) = -a/A + 2/i lalottj ^ (ai^-iGj + a2fi'2G^ + a^g^G^dx 



tanh ^a/ a + 2/i l^lgT^ y (ai/iiG{ + ^2^2^^ + a3/i3G^)(ix 
sinh(^VA + 2/i|a|o(T-t) 




an 



cosh I a/A + 2/i l^lgT 



(aiXiGj + a2-'^2G2 + 0:3X3^3) dudt, 
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cosh (-\/A + 2/i |q;|qT) 
X J (^aiUi (x, T) G{ + a2U2 (x, T) + 03^3 (a;, T) j dx, 

E2j{I){a) = -^/JI\a\Q j (^-\a\lgjG^j - aj{aigiG{ + a2g2Gi + asgsGl)^ dx 

— tanh (^/I |q;|qT) J (^— \a\^hjG-'j — aj[aihiG{ + a2h2G-'2 + a^h^Gi)^ dx 



T 

+ 




an 







- ^)) 


cosh (^y/JI 


\a\ 


\oT) 



X (^laloXjG^^ + aj{aiXiG{ + a2X2G{ + 03X3^0 j dudt, 



— Qfj ^aiUi (x, T) + 0:2^^2 {x, T) Gi + a-^u^ (x, T) j ^ (ix. 

Note that i?^^- and i?2j in Lemma [T] depend on u{x, T) instead of the data / = X, g, h). 
Therefore, in general these terms are unknown. However, our observation is that with |a| 
large, Elj and i?2j relatively small in comparison with Eij and E2j, and can be relaxed 
when computing the integrals J fjGdx which are the Fourier coefficients of f/s. So we 

introduce some convenient notations. 

Definition 1 (Information from data). For I = {ip, X, g, h) and a G such that a ■ a < 0, 
denote (using notations of LemmaU^. 

lo, tfD,{I){a).D2{I){a)=0. 
Definition 2 (Fourier coefficients). For a = (01,02,0:3) G C'^, w E L"^ {Q), denote 

J-'(w) (o) = / w{x)G{a, x)dx = / w(x) cos (oiXi) cos (02X2) cos (03X3) (ix. 



Note that any function w G L^(fi) admits the representation 

w{x) = fi:(m, n,p)J-' (w) (mvr, ?T,7r,p7r) cos(m7rxi) cos(n7rx2) cos(p7rx3), (5) 

m,,n,p>0 

where K{m,n,p) := (1 + l{„^o})(l + l{n5^o})(l + l^o})- 

As we explained above, we hope to approximate F{fj){a) by Hj{I){a) with |o| large. 
To do this, we need some lower bounds on \Di{I){a)\ and \D2{I){a)\. We shall require the 
following assumptions on and T. 



4 



(Wl) There exist A(v?) G (0,T) and C{ip) > such that either i^{t) > C{<^) for a.e 
t e (0, A(y?)) or ip{t) < -C{ip) for a.e t E (0, A(y?)). 



W2) T > 9 may J— _ ^ I or (W2') T > 12^5 x maxi— ,^^=l. 

Remark 1. If is continuous at t = then the condition (Wl) is equivalent to ^p{0) ^ 0. 
The conditions (W2) and (W2') mean that the observation time must be long enough. The 
condition (W2) is enough for the uniqueness, while the stronger condition (W2') is required 
in our regularization. These conditions should be compared to similar conditions (2.7) and 
(2.8) in 0. 

Theorem 1 (Uniqueness). Assume that (Wl) and (W2) hold true. Then the system (Hj^y 
has at most one solution 



f) e {{C\[0, T];L\Q)) n L\0, T; H\Q)))\ (L^Q))') . 



The main point in our regularization is to recover T{f){a) with |a| small from approx- 
imate values of J^{f){a) with |a| large. As in [7J we shall use the Lagrange interpolation 
polynomial. 

Definition 3 (Lagrange interpolation polynomial). Let A = {xi, X2, x„} be a set of n 

mutually distinct complex numbers and let w be a complex function. Then the Lagrange 
interpolation polynomial L [A; w] is 

n / 

Z- Xk 



l[am{^)=j:[ii^—v]^(-^)- 



Theorem 2 (Regularization). Assume that {u^, /°) is the exact solution to the system (Hj^ 
with the exact data Iq = {(p^ ^h^), where the conditions (Wl) and (W2') are satisfied. 
Let 6 > 0, consider the inexact data Is = {(f'^, X^, g"^, h^) such that, for all j G {1, 2, 3}, 



.,<eA\X'^-X^A\^,<eM-g^,\\^,<eM 



Construct the regularized solution /J from = {(f"^ , X'^ , g"^ , h^) by 

fj= K {m,n,p)J-'e {m,n,p) G {rmT,mi,p7i) 

<m, n, p < Ti; 

where 

V 60 ' 60 
Br, = {±(5r,+j):j = l,2,...,24rJ, 
J^e{fn,n,p) = L[Br,,Hj{Is) {-i.,mT,pTT)]{im7T) . 

Then we have, for all j G {1,2,3}, 
(i) (Convergence) G (R^) and /° in L'^{Q) as e 0. 
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(a) (L"^ -estimate) If fj G then fj — )■ fj in H^{Q) and there exist constants Eo > 

and Co > depending only on the exact data such that for all e G (0,£:o); 

\\f!-f!LHn,<Co{\n{e-')y\ 

(Hi) (H'^ -estimate) If f^ G H'^ {^) then for all e G {0,eo), 

\\fj-fj\\HHn)<Co{\n{e-')y\ 

Remark 2. In our construction, the convergence in H^{Q) is not expected even if fj G 
C^^iQ) since df^/dn = ondQ. 



3 Uniqueness 

In this section we shall prove Theorem [H We start with the proof of Lemma [1] by using the 
argument in [B]. 

Proof of LemmaUl Fix a = (01,02,0:3) G C"^ with \o\l = — (af + q;| + a|) > 0. For any 
j G {1,2,3}, getting the inner product (in L^(n)) of the k-th equation of the system ([1]) 
with Gl {k = 1,2,3), then using the integral by part and the boundary conditions and 
(jl]), we have 

72 



^ J u kG^dx + fi[ol + ol + oj) J UkG{dx 

XkGidio-^{t) I fkGidx. {A, 



an 



Multiplying {\Ak\j with OjOk, and then getting the sum of /c = 1, 2, 3, we have 
d^ 



dt^ 
o 



{ojOiUiGi + Oj02U2G2 + OjO^u^G-^) dx 

— (A + 2/i)|a|Q j {ajOiUiG{ + Oj02U2G^2~^ '^j'^'i'^'i^^'ijdx 
n 

J [ojOiXiGi + aja2X2G-!2 + aja-sX^Gi) doj 

+ '^(t) j {ajOifiGl + aja2f2Gi + aja-sfsG^) dx. (6) 
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Choosing k = j in ( [AfcD , then multiplying the result by |a;|o, and adding to ([H]), we obtain 
J (^\a\lujG + (^ajaiUiGl + aja2U2G{ + ajasUsGi^^dx 

= J ^lalgX^G^- + (^ajaiXiGl + aja2X2G{ + ajasX^Gi^^duj 
an 

+(f (t) J (^\a\lfjG^j + {ajaifiG{ + aja2f2Gi, + ajasfsG^^dx. (7) 
n 

We can consider (jH]) and ([7]) as the differential equations of the form 

y" (t) - riS (t) = hit), (8) 
where r] > is independent of t. Getting the inner product (in L^(0,T)) of ([8]) with 

T 

- y' (0) tanh (r^T) - (0) + ] y (T) = [ h (t) "'"'^^l^'f dt. (9) 

cosh [ril ) J cosh (r^i j 



Applying (E]) to (jH]) and ([7]) with rj = y/X + 2yU \a\Q and rj = y/]I\a\Q respectively, we get 

E,j (/) (a) + E*,^ (a) = -^D^ (/) (a) / a, + a2f2Gi + a^hG^) dx, 

— a n J 



'0 Q 



^2, (/) (a) + («) = -^D2il)ia)x 



— \a 







X / — a 



lo/jG^. - aj{aifiG{ + ^2/2^^ + asfsG^^dx. 



It follows from the latter equations that 

y A (/)(«) ^2 (/)(«) 

n 

This is the derised result. □ 

The following lemma gives a lower bound for \Dj (/) {a)\ (defined in Lemma [1]) when ip 
satisfies the condition (Wl). 

Lemma 2. Let ip E (0, T) satisfies the condition (Wl ) then there exists a constant R (ip) > 
such that for all a = (ai, 02, a^) G C'^, a ■ a < and \a\Q> R {p), 

|Z),(/)(a)|>lKC(^),„m{-L.^} /„r, = 1.2. 
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Proof. Denote ki = ^/X + 2fj, and /c2 = y/Ji- By the triangle inequahty, one finds that 

sinhl kj \a\n.(T — t) 



(/) (a) 



\a 



> 







fit) 



cosh (/cj lalgT) 
sinhl /cj |q;|q(T — t) 



dt 



A(V) 



cosh (fcj lalpT) 
cosh (/cj lalo?") 



dt 



dt 



We have, with |q;|o large, 



^ sinhf fc,- laLfT — t) 



A(^) 



cosh [kj \a\QT) 



dt 



< 



< 



sinh — A (1^9) 



cosh (/cj lalo^) 



4max{A;i, A;2}|«|o 
On the other hand, the condition (Wl) implies that 



^^'^'^ sinhf fc,- \a\AT — t) 



dt 



CM. 

kj |q;|q 



cosh {kj {qIqT) 
coshf kj I^Iq^T — A( 



cosh {kj IoIqT) 



> 



sinh^/cj |a|o(7' — i)^ 
cosh {kj \q.\qT) 

C{f) 



dt 



2kj\a\o 



with |a|o large. The desired result follows immediately from the above inequalities. 



□ 



The proof of the uniqueness below follows the argument in [7J. We shall need a useful 
result of entire functions (see e.g. [3j, Section 6.1). 

Lemma 3 (Beurling). Let (p be a non-constant entire function satisfying the condition: there 
exists a constant k > such that (r) < ke^ , for all r > 0, where (r) = {|0 (2;)! : \z\ = 
r}. Then 

In \(h(r)\ 
limsup '^^ ^' > -1. 

r— >oo 1" 

Proof of TheoremUl Suppose that {u^, f^) and {u^, f^) are two solutions to the system (PH]) 
with the same the data / = {ip,X,g, h). Then (u, /) := {u^ — m^, — /^) is a solution to 
(dMD with data ((^, 0, 0, 0). We shall show that (m, /) = 0. 

Assume that / 7^ 0, namely fj 7^ for some j G {1, 2, 3}. For any r;,,p G N U {0}, let us 
consider the entire function 



z I—)- Tpn^p {z) = I fj{x) COS {izxi) COS (^71X2) COS {pnx^) dx. 



Because 

rlih 

{imn) = / ixifj{x) sm{rmTXi) cos{nTCX2) cos{pTTX3) dx 



dz 

and {sin (mvrxi) cos (n7rx2) cos (p7rx3)}^gpj „pgpj^lQl is an orthogonal basis on there 
exist some (no,po) such that ipno,po is non-constant. 
On the other hand, recall from Lemma [1] that 

^.0.0 w = j mG{a..x)dx = :d7(7)(^+:d^(7h^' (10) 



where ar = (ir, no7r,po^)- Fix e > such that mm{Ty/ X + 2/i, T^/JI} > 2 + e. Then it 
is straightforward to see, from the explicit formulas in Lemma 1 and the lower bound in 
Lemma [21 that 

|^;-(«r)| < Cle-(l+^)^ l^e(/)K)| > Car, e = 1,2, 

with r > large, where Ci > and C2 > are independent of r. Therefore, it follows from 
(fTO|) that 

for r > large. This yields 

limsup < —(1 + e) 

r— >oo r 

which is a contradiction to Lemma [31 Thus / = 0. 

Now following the proof of Lemma [T] up to the equations and ([7D (a^ := al + al + 
needs not be negative at that time) we obtain 

d^ f 

+ (A + 2fi) c? J (^ajaiUiG{ + aja2U2G{ + ajasU^Gi^dx = 0, (11) 
n 

n 

+ fxa"^ J (^a'^UjG-'j — [ajaiUiG{ + ajaaMaG^ + aja^UsG-^^dx = (12) 
n 

for all a G C^. Note that y = is the unique solution to the differential equation 

y"{t)+r]^y{t) = 0, 
y{0)=0, 
y'{0) = 0, 

where 77 > independent of t. Applying this to (ITT]) and ( IT^ with a G we get 

l{a,a.u.Gi + a,a2U2Gi + a,a,u.GQdx^0 



dt^ 



and 



G — {oijaiUiGi + aja2U2G2 + aja^u^^G-i^ jdx = 0. 



Adding the latter equations we obtain 



(13) 



To get the same equation to (fT^ with a = 0, we can simply use identity with a = 
to get 

Uj{x, t)dx = 0. 



Since 



we have 





= = 




i=0 




J Uj (x, .)d 




n 



d 

dt 



Uj [x, t) dx 



t=Q 



(14) 



Putting ( lT3jl and f[T4j) together, we arrive at 

UjG = / .) cos(aia;i) cos(Q;2a;2) cos{asX3)dx = for all a G 



This implies that u = because {cos(m7ra;i) cos(m7rx2) cos{mTXs)}m,n,p>o forms an orthogo- 
nal basis of □ 



4 Regular izat ion 

Lemma [Hand Lemma [2] allow us compute J^{f){a) approximately with |a| large. To recover 
J-'{f){a) with \a\ small, we shall need the following interpolation inequality. This is an 
adaption of Lemma 4 in . The only difference is that we are working on a 3-dim problem, 
and hence we require more interpolation points. 

Lemma 4 (Interpolation inequality). Let Br = {±Zj\zj = 5r + j, j = 1, 2, . . . , 24r} for some 
integer r > 50. Let u : C ^ C be an entire function such that \u) {z)\ < Ae\^\ for some 
constant A. Then for any function g : Br ^ C one has 

sup \u (z) - L [Br, g] iz)\<A e'"' + 48 r e^^*^ sup \u {z) - g {z)\ . 

\z\<-Kr ZS-Br 

Proof. Fix 2; G C with \z\ < nr. We have the following residue formula at 49 simple poles 
{z}UBr, 



C={xeC:|x|=50r} 



cu [x) Y\ j 

X — Z x'^ — z'^ 

j=l J 



dx = 2711 (^u (z) — L [Br, uj] (z)^ . 
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This gives the simple bound 



50rsup|M^i4ff 



oj{z)-L [Br, u] {z) I < 50r sup <^ ^ ^ . (15) 





1 




- z'^\ 



Now we bound the right hand side of (1151) . For x G C, |x| = 50r, we have \uj {x)\ < Ae^^^, 
\x — z\ > (50 — tt) r and 



We shall show that 



\x 



z]-jj^i50rr-z^ 



24r / \ 2 9 

n Trr + 27 50 — TT ^i 



(vrr)^ + X 
(50r)^ — X 

[O, (29r)^]. Using Jensen's inequality we get 



Note that the function x ^-^ v (x) = In [ — ^ ) . v (x) is increasing and concave in 



j=l /i=l y=4(/i-l)r+l J h=l y j=4(/i-l)r+l 

= 4r ^ 1^ (2Ah + IG/i^ + y ^ + (3 + 4/i) ^ + ^ 
< 4r ^ i; M 24/i + IG/i^ + j + (3 + 4/i) — + 



6 X 50^ 



4r In 



vr2 + 24/. + 16/.2 + f + ^ + ^ 



< — 51r + In 



6x50^ 

50 -tt' 



50 

Thus ( 1T7|) holds. Replacing f lT^ and f lT7|) into ffT^ we obtain 

|a;(z) -L[5^,a;] (2)1 < Ae-^ ^z e C, \z\ < nr. {U 
Next, we observe that 

24r / 2 _ 2 \ , 

L [Br, u] (z) - L [Br, 9]iz)=Y.[U ) ^ (^^O - ^ (^^O) 



24r 



v2 ^2 
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Therefore 



24r 



L[Br,uj]{z)-L[Br,g]{z) < 2a^ 



n 



(19) 



where a = sup ja; (z) — g {z)\. 

Now we bound the the right hand side of (1191) . We have 



n 



< 



< 



n 



Zk 



Zj -\- Zj^j Z]^ \Zj Z}^\ 



n 



Zk 



Zj Zk 



< 



n 

(5r + 2) (5r + 3) . . . (29r) 
(j-l)!(24r-j)! 



(5r + 2) (5r + 3) ■ ■ ■ (29r) ^ 

(12r- l)!(12r)! -• "^l^J- 

Since J (1) < e^" and, by direct expansion, 

J(r + 1) (29r + 1) (29r + 2) . . . (29r + 29) 



Jir) 



< 



(5r + 2) . . . (5r + 6) ■ [(12r + 1) . . . (12r + 11)]^ ■ (12r) ■ (12r + 12) 
29^9 



< e 



30 



55 ■ 12^4 

we conclude that J (r) < e^"'" for all r > 1. Thus (fT9l) reduces to 



L[Br,tu] {z)-L[Br,g] (z) 



< ASrae 



30r 



The desired result follows from fllSp . (120]) and the triangle inequality. 



(20) 
□ 



The last of our preparation to prove the regularization result is the following useful lemma 
on the truncated Fourier series. The proof is elementary and we refer to Lemma 5 in ^ for 
details (in [S] the authors dealt with the 2-dimensional case, but the proof for 3-dimensional 
case is essentially the same). 

Lemma 5. For each w G L^(fi) and r > 0, we define 

Tr{w){x) = K{m,n,p)J^{w){m'7i,n7i,p7i) cos{m7TXi) 008(1271x2) cos^pnx^) 

0<m,n,p<r 

Then 

(i) Triw) w in L^iVt) as r —)■ 00. 
(a) If w & H^{fl) then Tr{w) w in H^{Q) and 



\w\ 



(ill) If we H\n) then 



|r^(w) - < ^ 



\w\ 
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We are now ready to prove the main result. 

Proof of Theoreml^ We shall first estimate the error /J — r^^/j*, where Tr is defined in 
Lemma Then we compare F^^ /° and /g by employing Lemma El The conclusion follows 
from the triangle inequality. In the following, C, Ci, C2 always stand for constants depending 
only on the exact data, and saying that e > small enough means that e G (0, eo) for some 
constant £0 > depending only on the exact data. 

Step 1: Bound on |F? (m, n,p) — J-" (/j*) (rnvr, n7r,p7r)| for m,n,p G [0,r^]. 
Applying Lemma |1] to g {z) = Hj {z, n7r,p7r) and 



U)[Z 



fj {xi, X2, X3) COS (—izxi) COS {nTcx2) COS (pvrxs) dxidx2dx3 



we find that for m,n,p & [0, r^] and e > small enough 

\Fj {m,n,p) — J" (/°) (mvr, n7r,p7r)| 



[im-K] 



< 



CO [irrni , 



l/iLnn)^"''^ +48^.^30'-^ sup \u{z)-g{z)\. 



(21) 



The error \u {z) — g {z)\ can be bound by direct computation. Indeed, for z G Br^, we 
have 



\uj{z)-g{z)\ 



= (fj) i-iz,mT,pTi) - {z,nTC,pTc)\ 
< {fj) i-iz,mT,pTT) - {z,mT,pTc)\ 
+ I Hj {z, nn, pn) — Hj {z, nn, pn) | 
Elj {—iz,mi,pTc) {—iz,mT,pTc) 



< 



+ 



Di (Jo) {-iz, mr, pir) D2 (/q) {-iz, nvr, pn) 



+ 
+ 



Eij {h 



-iz, nn 



pn) Eij{I^){-iz,nn,pn) 



Di (Jo) {-iz, nn,pn) Di {!,) {-iz, nn, pn) 
E2j (/o) {-iz, nn, pn) E2j {h) {-iz, nn, pn) 



E>2 {h) {-iz,nn,pn) D2 {h) {-iz,nn,pn) 
It is straightforward to see from the explicit formulas of Eij, E2j, Elj, E2j that 



\Eij {!,) {-iz,nn,pn) - Eij (Jq) {-iz,nn,pn)\ < Cir^^ e^^''' e, 
\E2j{Ie){-iz,nn,pn) - E2j{Io){-iz,nn,pn)\ < C^r^ e^^*"^ e. 



\Eij{Io){-iz,nn,pn)\ < CiV^e 



3 „29rr 



\E2j (/o) {-iz,nn,pn) 
\EIj {-iz,nn,pn) 



< CirU 



3 „29n 



< 



E;^{-tz,nn,pn)\ < -j^^- 



(22) 
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On the other hand, Lemma [2] ensures that C2T£ < \Dj (Jq) {—iz,mT,pTT)\ < Cirl and 

\Dj{Ie){-iz,n'R,p'K)\ > \Dj (Iq) {-iz,mT,pTT)\ 

- \Dj (4) {-iz, nvr, pn) - Dj (Jq) {-iz, rni, pn) \ 

for all z G Bj.^ and e > small enough. Therefore, from (I^T]) and fl22p we get 
|F/(m,n,p) - J^(/°) (m7r,n7r,p7r)| < \\f^\\^^^^^^e^ + 48 x 



1 / /=- ■ ^ „ _„\ 1 / /IT ™ _„\ / 



e 5 



Moreover, the condition (W2') gives 

i-<i(v/5VIT2^T-59), i<-L(ygVAr-59). 



60 60 V ^ " / ' 60 60 

Thus we obtain 

|F;(m,n,j9) - J^(/°) (m7r,n7r,p7r)| < c(ln (23) 
for all £ > small enough. 

Step 2: Conclusion. 

(i) Using the Parseval equality and the estimate fl25]) one gets 

W^l = K(m,n,p) |F/(m,ra,p) - (m7r,ra7r,p7r)|l 

< m, n, p < 

< C{\n{e-^)fh^o. (24) 

Since ||//-r,J°||^,(^) ^ by (121 and llr.J" - /"H^.^^^ ^ by Lemma we obtain 
||/| — /j'||L2(!^) — )■ as £ — !■ by the triangle inequality. 

(ii) Assume that /° G H^{Vl). Then by Lemma O one has HrrE/j — /?°ll_H'i(n) ~^ ^ ^^"^ 
[[Fr^/j' — /j ||j;^2(j7) ^ C(ln(£:^^))^5. Combining the latter estimate with fl24|) . we conclude 
using the triangle inequality that 

\\r-f'\Un)<c{}^{e-')y" 

for e > small enough. 

Moreover, using (l23i) and the Parseval equality for iiT^ we have 



P0||2 

= K (m, ri,p) (1 + 7r^(m^ + + p^)) |FJ (m, n,p) — (/°) (mTT, ?7,7r,p7r) 1^. 

0<m,n,p < 

< C{ln{e-^)Yh^^ 0. (25) 
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Thus — llji^i^Q) by the triangle inequahty (in if^-norm). 
(iii) Assume that fj G H'^{Q). Then by Lemma [5] one has 

This estimate together with ( 125|) and the triangle inequality yield 

for e > small enough. □ 



5 Numerical example 

In this section, we test our regularization process in an explicit example. Choose fi = 1, 
X = —1, T = 30 and consider the system (PS]) with the exact data Iq = (^9°, g^, h^) given 

by 

(^0 {t) =237r2 cos (vrt) , /i^ = /i^ = = 0, 
X'^ (xi, X2, Xs, t) = COS (vrt) ^— 47r sin (27ra;2) sin (27ra;3) rii 

— 271 sin (47rxi) sin (27rx3) n2 — 27r sin (47rxi) sin (271x2) , 
X2 {xi,X2, X3, t) = COS (7rt) ^— 27r sin (47rx2) sin (271x3) ^1 

— 4:71 sin (27rxi) sin (27rx3) n2 — 27r sin (27rxi) sin (47rx2) ^3 j , 
X3 (xi, X2, X3, t) = cos (nt) (^—27!- sin (27rx2) sin (47rx3) rii 

— 271 sin (27rxi) sin (471x3) ^2 — 47r sin (27rxi) sin (27rx2) n^j , 

(xi, X2, X3) = sin (47rxi) sin (27rx2) sin (27rx3) , 
g2 (xi, X2, X3) = sin (27rxi) sin (47rx2) sin (27rx3) , 
5^3 (xi, X2, X3) = sin (27rxi) sin (271x2) sin (471x3) ? 

Note that in this example the conditions (Wl) and (W2') are satisfied. The exact solution 
of the system (PHj) can be computed explicitly 

Ui (a^i, X2, X3, t) = cos (nt) sin (47rxi) sin (277x2) sin (277x3) ? 

U2 (xi, X2, X3, t) = cos (7rt) sin (27rxi) sin (47rx2) sin (27rx3) , 

■U3 (xi, X2, X3, t) = cos (7rt) sin (27rxi) sin (277x2) sin (477x3) , 

fi (xi, X2, X3) = sin (47rxi) sin (277x2) sin (277x3) ? 

/a (xi, X2, X3) = sin (277Xi) sin (477x2) sin (277x3) , 

/g (xi, X2, X3) = sin (277Xi) sin (277x2) sin (477x3) • 
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Now, we consider the disturbed data, for n G N, j G {1, 2, 3}, 

— StT cos [TTt) / 

X" = X° + ^ sin (2n7rx2) sin (2^71x3) rir 



+ sin {2mTXi) sin (2r;,7rx3) n2 + sin {2mTXi) sin (2r;,7ra;2) na 



9- =9' + 



sin (2n7ra;i) sin (2?7,7ra;2) sin {27111x3) 



3 ' 3 



The disturbed solution of the system (PS]) with the disturbed data is 
^ Q ^ cos (vrt) sin (2n7ra;i) sin {2mTX2) sin (2n7rx3) 



V ^ 



Q (— 1 + 12n^) sin (2n7rxi) sin (2n7rx2) sin {2mixs) 



23Vn3 

We can see that the error of the data is small 



(in fact, they even converge in the uniform norm). However, the error of the solution is large 
since 

,011 V2 /I44n4 - 24n2 + 1 



V! fjWL^n) 4 Y 529n3 ^ °° 

as 77, — )■ cxD. Thus the problem is ill-posed and a regularization is necessary. 

Now, we apply our regularization procedure for e = 0.01 and the disturbed data with 
n = 10. The resulting regularized solution is 

fl (xi, X2, X3) ~ 0.035 — 0.063 cos (irxi) — 0.156 cos {nx2) — 0.455 cos (7ra;3) 
+ 0.033 cos (vrxi) cos (7rx2) + 0.027 cos (vrxi) cos (ttx^) 
+ 0.15 cos (vrx2) cos (vrxs) — 0.005 cos (ttxi) cos (vrx2) cos {1TX3) . 

The error between the regularized solution and the exact solution is 

l|/i-/?llla(.)- 0-273. 
To see the effect of our regularization, note that the corresponding disturbed solution (with 

1 1 2 

data error e = 0.01) causes an extremely large error \\fj^ — fj\\j^2(^Q-^ ~ 3.4 x 10^. 

Acknowledgments. The work was done when P.T.Thuc was a Master student in Ho 
Chi Minh City University of Education. 
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